


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1977 


Numerical modelling of the evolution of 
monsoon circulation along 80 degrees E. 


Schwenke, George William 


Monterey, California. Naval Postgraduate School 


http://ndl.handle.net/10945/18218 


Downloaded from NPS Archive: Calhoun 


: Calhoun is the Naval Postgraduate School's public access digital repository for 
/ (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist : Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

; | LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


NUMERICAL MODELLING OF THE EVOLUTION 
OF MONSOON CIRCULATION ALONG 80'E 


George William Schwenke 





NAVAL POSTGRADUATE SCHOOL 


Monterey, California 





THESIS 


NUMERICAL MODELLING OF THE EVOLUTION 
OF MONSOON CIRCULATION ALONG 80°E 


by 


George William Schwenke 


March 1977 


Coasb Chang 
Thesis Advisors R. T. Williams 





Approved for public release; distribution unlimited. 


T178051 





a a rr SAP a ee 
SECURITY CLASSIFICATION OF THIS PAGE (When Date Entered) 


REPORT DOCUMENTATION PAGE 


— 


4. TITLE (and Subtitie) 


Numerical Modelling of the Evolution 
of Monsoon Circulation Along 80°E 






READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


{ 3. RECIPIENT'S CATALOG NUMBER 


3. TYPE OF REPORT & PERIOO COVERED 


Master's Thesis 
977 


REPORT NUMBER 















Mie 
6. PERFORMING ORG. 


7, AUTHOR(e) 6. CONTRACT OR GRANT NUMBER(@) 


George William Schwenke 


10. PROGRAM ELEMENT, PROJECT, TASK 
AREA & WORK UNIT NUMBERS 


12. REPORT OATE 
March 1977 


13. NUMBER OF PAGES 


40 


SECURITY CLASS. (of thie réport) 








9. PERFORMING ORGANIZATION NAME ANO AODORESS 


Naval Postgraduate School 
Monterey, California 93940 : 


















CONTROLLING OFFICE NAME ANO ACORESS 


Naval Postgraduate School 
Monterey, California 93940 























MONITORING AGENCY NAME @ ADDRESS(If different from Contraiting Oflice) 


Naval Postgraduate School 
Monterey, California 93940 





Uitheeeais salmaed 


DECL ASSI FIC ATION/ DOWNGRADING 
SCHEDULE 





Se. 





'§. DISTRIBUTION STATEMENT (of thte Report) 


Peemoved=tOr public release; distribution unlimited. 


a 


17. DISTRIBUTION STATEMENT (of the abatract entered in Block 20, If different from Report) 


18. SUPPLEMENTARY NOTES 


19. KEY WORDS (Continue on reveree aide if neceeeary and identify by block number) 


Monsoon General circulation model 
Zonally-symmetric 


Parameterization 
Easterly Jet 





20. ABSTRACT (Continue on reveree side if neceeeary and identify by block number) 


The UCLA general circulation model developed by Arakawa and 
Mintz (1975) was used to study the evolution of the northern 
Summer monsoon during a 180 day period commencing in mid January. 
The global model, modified to a five degree latitude by five 
degree longitude grid resolution, was run for five vertical levels| 
in a zonally svmmetric mode along 80°E to permit the extended 
computational period. Thermal forcing was simulated by 










DD , Pee 1473 EDITION OF 1 NOV 68 1S OBSOLETE 


(Page 1) 1 SECURITY CLASSIFICATION OF THIS PAGE (When Dare Entered) 





yp a SS 
Foc uMty CL aseiric ation OF TemlS PAaGE(Wren Neta Entared. 


continuously adjusting computed temperatures at each level toward 
the radiative equilibrium temperature for that level as esti- 
mated by Palmén and Newton (1967). These equilibrium temper- 
atures were also adjusted during the 180 day period from the 
January mean to the July mean. The evolution and location of 
the major zonal components of the monsoon were reasonably simu- 
lated. The results indicate that a zonally symmetric model is 
suitable for studies on the planetary-scale features of the 
monsoon with the realization that exclusion of zonally-asymmetri 
phenomena from such a model limits its accuracy and usefulness 
in considering the effects of such events on the monsoon system. 


DD ekg n 1473 
Es * an i 
S/N 0102-014-6601 2 SECURITY CLASSIFICATION OF THIS PAGECMREN Date Entered) 





Approved for public release; distribution unlimited 


Numerical Modelling of the Evolution 


of Monsoon Circulation Along 80°E 


by 


George William Schwenke 
Lieutenant, United States Navy 
B.S., University of Wisconsin, 1971 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCTENCE IN METEOROLOGY 


“.from the 


NAVAL POSTGRADUATE SCHOOL 
March 1977 





TABLE OF CONTENTS 


ir. INTRODUCTION - - - = = = = = - = 2 = a we e = 


fl ewbscee te itehP OF THE NUMERICAL MODEL 


A 
B 
. 
D 
E 


3 
i] 
3 
i] 
5 
| 


BASIC MODEL - - - - - -« - - - = «© = = - = = 
THE GRID SYSTEM - - - - - - - = - - - - = - 


PARAMETERIZATION OF SUB GRID-SCALE PROPERTIES 


PiARMTe Senne =e 
INTEGRATION SCHEME AND INITIAL CONDITIONS - 


mee RESULTS - - = = = = = 2 se ee ee ee eee 


IV. CONCLUSIONS - - - - - - - - = = = = - - = - - - 


LIST OF REFERENCES - - - - - - - = - - - -+ = = = - = 


INITIAL DISTRIBUTION LIST 


f 
f 
j 
a 
J 
¢ 
i 
f 
’ 
g 
f 
J 
f 


LZ 
13 
14 
14 
1) 
16 
1 
58 
a0 





ILE 
Zn 


San 


Zonally-averaged temperatures for January 


Vertical! Vemels - 


LIST OF FIGURES 


= 


Horizontal distribution of variables 


Zonally-averaged temperatures for July - 


Observed zonal geostrophic wind along 80°E 


Zonal 
Zonal 
Zonal 
Zonal 
Zonal 
Zonal 
Zonal 
Zonal 
Zonal 
Zonal 
Zonal 
Zonal 
Zonal 


Zonal 


WW Component 
u component 
u component 
u component 
u component 
u component 
u component 
Gemperacure 
temperature 
Bemperature 
tenpcracure 
Eenperacure 
eempera cure 


Genperacture 


for January 


for February 


for 
for 
LOT 
for 
satel g 
for 
for 
for 
for 
On 
LO % 


ror 


March - 
hou. 
May = = 
ete a 
ay = 
January 
February 
March - 
April - 
May - - 
Jie Ss 


Sly =a 


13 
7a 


ed 
‘ss 
ZS 
24 
es 
26 
oy 
Ze 
oo 
30 
31 
G4 
55 
34 
ope 
36 
o/ 





eh KHOX LeRapy 
POSTGRADUATE SCHon: 


ABSTRACT 


The UCLA general circulation model developed by Arakawa 
and Mintz (1975) was used to study the evolution of the north- 
ern summer monsoon during a 180 day period commencing in mid 
January. The global model, modified to a. five degree lati- 
tude by five feeres longitude grid resolution, was run for 
five vertical levels in a zonally symmetric mode along 80°E 
to permit the extended computational period. Thermal forc- 
ing was simulated by continuously adjusting computed temper- 
atures at each level toward the radiative equilibrium tem- 
perature for that level as estimated by Palmén and Newton 
(1967). These equilibrium temperatures were also adjusted 
during the 180 day period from the Jemwesay mean to the July 
mean. The evolution and location of the major zonal compo- 
nents of the monsoon were reasonably simulated. The results 
indicate that a zonally symmetric model is suitable for 
Studies on the planetary-scale features of the monsoon with 
the realization that exclusion of zonally-asymmetric phenom- 
ena from such a model limits its accuracy and usefulness in 


considering the effects of such events on the monsoon system. 
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I. INTRODUCTION 


The monsoon circulation has long been the subject of con- 
jecture and analysis by meteorologists because of its plane- 
tary scale and its environmental effects on almost half the 
world's population. While for many people the word "monsoon" 
is synonymous with the occurrence of heavy seasonal rainfall 
over India and portions of southeast Asia, this event is but 
one component of the total monsoon, which extends its influ- 
ence from the surface to the stratosphere over one third of 
the globe. 

Krishnamurti and Bhalme (1976) defined the broad scale 
monsoon system in terms of several parameters which include: 
(1) the monsoon trough, (2) the Mascarene sub-tropical high, 
(3) the low-level cross equatorial jet, (4) the Tibetan 
high, (5) the tropical upper-tropospheric easterly jet and 
(6) the monsoon rainfall. The principal time-mean circula- 
tion features of the northern hemisphere summer monsoon have 
been known and understood for some time, but detailed syn- 
optic analysis and accurate numerical simulation are fre- 
quently hindered by the lack of meteorological data, espe- 
Cially for upper air levels. 

The monsoon develops over southeast Asia as a whole, 
extending roughly from 60°E to 120°E. The reversal of flow 
pattern between seasons is most intense and persistent over 


India and extreme southeast Asia. These regions plus the 





Tibetan plateau in southwestern China may be viewed as the 
most active areas for the summer monsoon. The two components 
which affect upper levels in these areas most directly are 
the Tibetan high and the Indian jet. The Tibetan high is a 
warm core high whose formation is thought to be related to 
the thermal and mechanical effects of the Tibetan highland. 
Latent heat release due to shallow and deep cumulus convec- 
tion appears to play an important role in the maintenance of 
the upper anticyclone. 

The easterly upper-tropospheric jet near 150 mb is another 
interesting feature of the monsoon system. This jet has winds 
of roughly 80 to 100 kts with the strongest winds being found 
[omc o the west or the southern tip of India over the Ara- 
bian Sea. The jet forms in the monsoon month of June and 1s 
usually present until September. The proper simulation of 
these two related components of the monsoon system by a 
given numerical model is frequently considered to be a pre- 
requisite for further application of the model to monsoon 
Studies. 

The observed circulation systems and the theories proposed 
for their existence have Poneneced numerous synoptic and 
numerical studies of the monsoon circulation. The numerical 
models have ranged from simple barotropic regional models to 
complex global general circulation models (GCM). 

Holton and Colton (1972) used a simple linear one-level 
barotropic vorticity equation with the observed 200-mb sea- 


sonal mean divergence as the forcing to successfully simulate 





the seasonal mean flow pattern. They found, however, that 

an extraordinarily large friction was needed to damp the 
vorticity generation and keep the field in phase with the 
forcing function. Bellis (1975) and Pentimontti (1976) used 
a multi-level model formulated by Monaco and Williams (1975). 
This model was essentially similar to the UCLA general cir- 
culation model in use at that time except that parameteriza- 
tion of sub grid-scale processes and Separate calculations 
for boundary layer terms were not included. The model fairly 
accurately simulated the observed monsoon features except 
that phase shift problems were also noted. Hahn and Manabe 
(1973) and Washington and Daggupath (1975) have demonstrated 
the ability of the complex GCM, which includes parameterized 
heating and topography, to successfully simulate the northern 
Summer monsoon. These complex GCM's have some distinct dis- 
advantages in their application to monsoon studies because 
the large geographical domains monsoon simulations require, 
coupled with the complexity of the physics in these models, 
Peothe in Vast Computational time and computer core require- 
ments. 

In an attempt to incorporate the improved vertical reso- 
lution and surface terrain representation found in the GCM 
into a more tractable model for monsoon studies, Murakami 
(1969) noted that the main characteristics of the monsoon 
circulation are practically the same along a given latitude 
over the area between roughly 60°E to 120°E. The geography 


of this area is relatively uniform in that it consists of a 
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land mass to the north and oceans to the south. He thus con- 
cluded that one could introduce an idealized numerical model 
using a two-dimensional vertical plane along 80°E from the~ 
equator to the north pole. This model thus considers only 
the meridional circulation and its associated zonal current 
and allows no disturbances to be superimposed upon the zonal 
motion. Murakami developed a nine level zonally-symmetric 
model which included as dependent variables horizontal and 
vertical velocity components, pressure, temperature and a 
moisture variable. Sea surface temperature was prescribed 
south of 10°N with land to the north where smoothed topo- 
graphy was included. TWO factors which were disregarded in 
this model were the seasonal variation of the monsoon cir- 
culation and the interaction between the northern and south- 
ern hemispheres. The model was started from an atmosphere 

at rest and reached a statistical steady-state after 80 days. 

Murakami found that a zonally-symmetric model could accu- 
rately simulate many of the well known broad-scale features 
of the monsoon, including (1) low-level southwesterly flow 
between the equator and about 30°N, (2) the upper-level anti- 
cyclone and (3) a pronounced easterly jet near 150 mb on ie 
wiencinmmupoundary of the upper-level anticyclone. 

Pemiomune Intent @f this study to evaluate the effective- 
ness and accuracy of the UCLA GCM as developed by Arakawa and 
Mintz (1975) in a zonally-symmetric mode along 80°E. This 
Study will differ from the model developed by Murakami in 


that the UCLA GCM incorporates parameterized heating and 


ot 








cumulus convection and a pole-to-pole geographical domain. 
Additionally, improvements in computer systems have increased 
computational speed to the extent that it is possible to 
initialize the model to the observed January zonally-averaged 
temperatures with the atmosphere at rest, let the model attain 
Steady-state and then Mite ceate in time for 180 days while 
Simultaneously adjusting the zonally-averaged equilibrium 


temperatures from January to July. 


II. DESCRIPTION OF THE NUMERICAL MODEL 


A. BASIC MODEL 


The UCLA general circulation model, developed by Arakawa 
and Mintz (1975) and modified by Arakawa and Lamb (1976) was 
used as the basic model. The primary prognostic variables 
are horizontal velocity, temperature and surface pressure, 
governed, respectively, by the horizontal momentum equation, 
the thermodynamic energy equation and the surface pressure 
tendency equation. In this study the model was run for a 
vertical plane extending from 90°N to 90°S along 80°E with 
five poememmiicevyels as the vertical coordinates in the spher- 
ical coordinate model. Cyclic continuity was maintained on 
either side of the primary 80°E vertical plane, with compu- 
tations at the poles stabilized by the zonal symmetry of the 


model. 


daz 





Ee ihe eweRiID SYSTEM 


In the vertical the model was run with five sigma levels 
Se Gee, 15, and .15) interpolated to constant pres- 
Perewsourraces of 900, 700, 500, 300 and 150 mb. The sigma 
coordinate is defined as 

a? 


¢ = (2.1) 


where p is the pressure of the sigma level, P, is the con- 
stant tropopause height which was set to 100 mb, and 7m is 


the "terrain pressure" defined as 
T = Den Or) 


The surface pressure, Pe» iSescteroehoO0 mb initially. 


Prom (2.1) 2t tollows that 


o = 0 at p 


I 
a 
ct 


o=latp#= P. 


which are the vertical boundaries of the model. The bound- 
ary conditions at o = 1 and o = O are both oe = 0. 

The horizontal distribution of variables is illustrated 
in Figure 2. The grid points are spaced every five degrees 
in latitude. Values of Vey eabewmlatcdenor 1mplicit pomnts to 
the north and south of the (i,j) point, are averaged to 
mean a value for the (1,j]) point. The implicit u values 


are equal due to the zonal symmetry of the model as used in 


this study. 
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C. PARAMETERIZATION OF SUB GRID-SCALE PROPERTIES 


A model of the interaction of a cumulus cloud ensemble 
with the large-scale environment, as developed by Arakawa 
and Schubert (1974), is applied in the UCLA GCM used in this 
Study. 

In the model, the large scale environment is divided into 
the subcloud layer and the region above. The time changes of 
the environment are governed by the heat and moisture budget 
equations for the subcloud mixed layer and the region above, 
and by a prognostic equation for the depth of the mixed layer. 
In the environment above the mixed layer, cumulus convection 
affects the temperature and moisture fields through cumulus- 
induced subsidence and the detrainment from the cumulus 
clouds of saturated air containing liquid water, which evap- 
Oorates into the environment. In the subcloud mixed layer, 
GUumulus Gonvection does not act directly upon the temperature 
and moisture fields, but it affects the depth of the mixed 


layer through the cumulus-induced subsidence. 


D. DIABATIC HEATING “ 


The UCLA model is designed to use a detailed solar radia- 
tion balance routine to calculate the net diabatic heating. 
This radiation subroutine calculates incoming short wave 
radiation as a function of solar zenith angle and computes 
eller actors as reflection from clouds, re-radiation from the 


surface and the effects of moisture on this long wave 
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radiation. However, this routine is very time consuming and 
thus was not used in this 180 day study. In this zonally- 
Symmetric application of the model, the forcing function is 
the zonally symmetric observed thermal gradient, or equili- 
brium temperature T*. These temperatures, taken from the 
observations given by Palmén and Newton (1967) are used in 


the following formula to calculate the heating component: 


a5 


On= = Paes) | 


In this adjustment formula T is the computed temperature for 
the given latitude and sigma level; T* the equilibrium tem- 
perature for this same point, and t a specified adjustment 
time taken as two days in this study. The T* field is ad- 
justed each hour during the 180 day computational period to 
reflect the seasonal temperature change from January to July 
Witle maintaining the north-south temperature gradient. Fig- 
ures 3a and 3b detail the T* fields for January and July, 
respectively. A linear interpolation is used to determine 
the hourly adjustment applied.to the January field during the 
180 day period to increment it to the July values at the end 


Of the period. 


E. INTEGRATION SCHEME AND INITIAL CONDITIONS 


The time integration is comprised of continuous sections 
of 1 Matsuno (Euler Backward) step and 4 leap frog steps. 


The model was initialized to the January zonally-symmetric 


Jess 





Saullubrium temperatures (T*) and specific humidity. Initial 
sea surface and surface temperatures for 80°E were estimated 
from the global data prepared by Schultz and Gates (1972) for 
January. The model was run with a 5°-resolved smoothed topo- 
graphy derived from the actual topography along 80°E. The 
model was integrated for 60 days to bring it into equilibrium 
with the January T* fields, and was then run for the 180-day 
period of study while the T* fields were changed to simulate 
the seasonal radiation changes. All fields were interpolated 
to five constant pressure surfaces of 150 mb, 300 mb, 500 mb, 


790 mb and 900 mb. 


itl RESUS 


The evolution and location of the major zonal components 
of the monsoon were reasonably simulated by the model in this 
See uchys. 

Figure 5a shows the calculated u component field for 
January. This field is the resultant of the 60 day integra- 
tion during which the model reached a steady wind state with 
the January equilibrium temperatures. The calculated field 


ae SS) ay cua 


shows a maximum of approximately +18 m sec 
the northern hemisphere mid-latitudes, with a broad, weaker 
area of upper westerlies in the southern hemisphere. A maxi- 
mum in the easterly wind is noted in the tropical southern 
hemisphere, reflecting the trade winds. The January temper- 


ature field, shown in Figure 6a, closely resembles the Janu- 


ary equilibrium temperature, as expected, except that the 
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calculated temperatures in the tropics appear to be slightly 
cooler than the equilibrium temperatures. 

Uivemie the Mmext thirty days, from day 1 to day 30, the 
gradual change in equilibrium temperatures was begun. By day 
30, the u component field (Figure 5b) showed a slight decrease 
in the intensity of the northern hemisphere upper westerlies, 
with little change elsewhere in the hemisphere. The temper- 
meume field for day 30, shown in Figure 6b, indicates a 
fairly significant lower-level warming in the northern middle 
and high latitudes, with only slight cooling in the southern 
hemisphere. 

iremiexemsimty days Of integration, from day 30 to day 90 
(mid February to mid April), showed the onset of a signifi- 
cant positive u component near 10°N latitude. Figures 5c 
and 5d detail these fields. The v components in this area 
were also observed to be positive, implying a southwest low- 
level flow. Additionally, around day 75, an area of upper- 
tropospheric easterly winds is noted at about 10°N from 200 mb 
to 400 mb. This feature appeared during the ten day interval 
(the interval for which data was saved during the run) between 
day 70 and day 80. The northern upper level westerlies con- 
tinued to weaken slightly during this sixty day period, and 
the temperature fields (Figures 6c and 6d) indicate continued 
northern hemisphere lower level warming, with the 900 mb 
freezing level extending up to latitude 63°N by day 90 (mid 
mpril). 

The period from mid April (day 90) to mid June (day 150) 


shows a continued intensification in both the low level south 
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westerly flow near 15°N and the upper level easterly flow 
near 10°N from approximately 150 mb to 500 mb. A bi-modal 
lower level easterly wind maximum is noted in the southern 
Mmemisphere tropics near 10°S and 26°S by mid June. This 
feature first appears in April as a single maximum near 28°S, 
with the two maxima developing in May. The southern hemi- 
sphere westerlies intensified and moved northward to about 
30°S by day 150 (June) with maximum amplitudes on the order 
of 19 m sec}, 

The temperature fields for days 90 to 150, as shown in 
Figures 6d through 6f, indicate continued northern hemisphere 
warming and southern hemisphere cooling. 

The last thirty days of integration, from day 150 to day 
180, represent the northern summer monsoon circulation as 
Simulated by the zonally symmetric model. The southwest low- 
level flow is quite well developed, with a maximum amplitude 
eeel9) m sec } near 14°N. The easterly jet is also evident, 
with a maximum value of about 10 m sae * Neanels Neat  ls0tmb. 
The northern hemisphere westerlies are displaced to near 55°N, 
with values of 10 m sec} at 350 mb, while southern hemisphere 
-westerlies in excess of 20 m sec} are located near 300 mb 
eae. 2 Sasos 

Temperatures for day 180, as shown in Figure 6g, are in 
excess of 20°C at 900 mb for northern hemisphere land areas 
south of 54 N. These are warmer by 3 to 7°C than the equili- 
brium temperatures for these areas, and this is perhaps a 
result of the cumulus parameterization affecting the calcu- 
lated temperatures by warming the lower andmiddle troposphere. 
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ive CONGEUS TONS 


The zonally symmetric model appears to be suitable for a 
wide range of studies on the monsoon and its evolution in 
memes In this study, the major zonally symmetric fontuees. 
of the monsoon, such as the tropical upper-tropospheric 
easterly jet and the low-level southwesterly flow regime, 
were accurately simulated in position. The magnitudes of 
these features were not so accurately simulated, however, as 
the upper level easterlies and westerlies tended to be under- 
predicted and the low-level flow appeared to be over-predicted. 
Figure 4 shows the observed magnitudes of these features dur- 
ing July 1967; the easterly jet has winds in excess of 
30 m sec while this zonally symmetric simulation produced 
winds of only about 10 m soe > This problem appears to be 
related to the thermal forcing used in this simulation. The 
zonally-averaged equilibrium temperatures, used in this study 
as the mean thermal gradient and the forcing function for 
mayor zOmal Circulation systems, are not representative of 
the more intense gradients which exist over Asia. As such, 
the use of these T* fields should and did result in weaker 
gradients over the monsoon area and thus weaker zonal circu- 
lation patterns in both easterly and westerly flow. 

The excessive magnitude of the low-level southwesterly 
flow is not so easily understood. The observed magnitude is 


approximately 10 m sect; while the zonally symmetric model 
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indicates values of nearly 20 m sec’t. The discrepancy may 
De aeresult of the high topography to the north and its 
emrects Om the calculation of lower level thermal gradients, 
as well as its effects on cumulus parameterization and plan- 
etary boundary layer parameters. 
Future work with zonally symmetric models might center 

around one of the following changes to this study: 

(1) Increase the vertical resolution of the model 
Prom the present five tropospheric levels to nine or more. 

(2) Use an equilibrium temperature (T*) field which 
more accurately reflects the true gradients along the longi- 
meade under studys 

(oeeeseetne actual radiation subroutine as the ther- 
mal forcing rather than equilibrium temperature fields. 

(4) Evaluate the effects of southern hemisphere sea 


surface temperatures on the northern summer monsoon. 
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Figure 1. Vertical levels: Sigma (0) is the vertical coordinate. 


Values for horizontal velocity components, geopoten- 
tial and temperature are carried at the five sigma 
levels (dashed lines). Values for terrain pressure 


(tT) and vertical motion ao are calculated on solid 
line pressure surfaces. \ dt 
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Figure 2. The horizontal distribution 
of dependent variables around 
a oivens(?,9)) Grid point. 
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The time-and-longitude averaged temperature 
(T) in January as estimated by Palmén and 
Newton (1969). Values are in degrees K. 





The time-and-longitude averaged temperature 
(T) in July as estimated by Palmén and Newton 


(1969). Values are in degrees K. 
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Zonal u component for day 30 (February) 
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Figure 5c. Zonal u component for day 60 (March) 
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Zonal u component for day 90 (April). 
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Figure 5e. Zonal u component for day 120 (May) 
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Figure 5f£. Zonal u component for day 150 (June) 
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Figure 5g. Zonal u component for day 180 (July) 
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Figure 6a. Zonal temperature for day 1 (January) 
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Figure 6b. Zonal temperature for day 30 (February) 
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Figure 6c. Zonal temperature for day 60 (March) 
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Figure 6d. Zonal temperature for day 90 (April) 
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Figure 6e. Zonal temperature for day 120 (May) 
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Figure 6f. Zonal temperature for day 150 (June) 
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Figure 6g. Zonal temperature for day 180 (July) 
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